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We calculate adiabatic charge and spin pumping through a serial double quantum dot with strong 
Coulomb interaction, coupled to normal metal or ferromagnetic contacts. We use a real-time dia- 
grammatic approach in the regime of weak coupling to the reservoirs. In the case of weak interdot 
tunnel coupling we investigate the influence of tunnel-induced renormalization effects due to charge 
fluctuations on the pumped charge and spin. We show that tunneling through thermally excited 
states can play an important role in the strong interdot coupling regime. In particular, for ferro- 
magnetic contacts, both effects enable the generation of pure spin currents. Furthermore they can 
lead to an inverted spin-valve effect or even to the inversion of the transport direction going along 
with a diverging tunneling magnetoresistance. 

PACS numbers: 72.25.-b, 73.23.Hk 
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I. INTRODUCTION 



Pumping through mesoscopic devices is realized in the 
absence of an external bias by the periodic variation of 
certain system parameters in time, thus creating directed 
transport of electrons from one contact to the other. 
When this variation is slow compared to the character- 
istic lifetime of the electronic states on the device, one 
speaks of adiabatic punrpingi^r— In recent years, experi- 
ments on quantized pumping were performed, aiming at 
the realization of a quantum standard for the current 
or for the repeated initialization of coherent quantum 
states, useful for quantum operations.— ~— Fewer experi- 
mental realizations dealt with the role of quantum in- 
terference effects giving rise to a directed current due to 
time-dependent periodic fields . 10 i 11 Furthermore, when 
the system through which the charge is pumped is small, 
as say, quantum dots, the Coulomb interaction cannot 
be neglected because of the small capacitance of the sys- 
tem. This impact of electron-electron interactions on adi- 
abatic pumping through quantum dots, quantum wires 
and metallic islands has lately been considered in the- 
ory— _ — Particularly, charge pumping can be uniquely 
due to Coulomb interactions in the non-linear regime 2 ^ 
or due to interaction-induced renormalization effects^ 
Adiabatic pumping requires the time-dependent modu- 
lation of at least two of the system's parameters. In the 
work presented here we are interested in charge and spin 
pumping through a serially coupled double dot, where 
the energy levels of both quantum dots can be indepen- 
dently varied by time-dependent gate voltages. This was 
so far examined for noninteracting system o 24 ' 25 with nor- 
mal conducting leads, and further studies were performed 
in the high-frequency regimei 2 ^— The study of double 
dots is particularly interesting due to their complex in- 
ternal spectrum: the eigenstates of the double dot are 
differently coupled to the leads, and these couplings are 
effectively energy-level dependent. This results in new ef- 
fects for pumping, as demonstrated in this paper, which 
are controllable through the mere modulation of the bare 



dot energy spectrum. Importantly here, this choice of 
two independently tunable parameters, i.e., the spectra 
of the two dots, has an easier experimental access than 
any available choice of parameters in a single dot. The 
fact that adiabatic pumping acts as a useful spectroscopy 
toolj22 revealing system specific features which are not 
accessible through a static measurement, establishes one 
of our motivations to study adiabatic pumping through 
interacting double dots in various regimes and setups. 
Previous studies on transport through a double-dot sys- 
tem with Coulomb interaction in the static regime re- 
vealed charge fluctuation effects, related to the inter- 
play of the double dot's internal structure and the reser- 
voirs ; 29 ' 30 and double dots were studied in an interferom- 
eter setup; 3 - 1 - to name a few examples. Lately double dots 
have been extensively proposed and used as measurement 
devices: the read-out of spin properties is discussed in 
Ref. [32| and references therein, the Pauli spin blockade 
was used to investigate relaxation times 3 - 3 - where the in- 
fluence of spin-orbit coupling and the role of nuclear spins 
have been treated. Furthermore, a double quantum dot 
device can act as a noise detector^ In Refs. [35| transi- 
tions in the double-dot spectrum due to a bosonic envi- 
ronment were studied. 

In addition to the charge degree of freedom also the spin 
degree of freedom of the electrons plays an important 
role in transport, in particular, in the context of spin- 
tronics 3 ^ and applications in quantum computation^ 7 - 
The spin-valve effector— in systems containing differ- 
ently polarized fcrromagncts and its tunability, e.g., by 
a gate control 4 ^ are therefore of interest. In the study of 
quantum dot spin valves^ 2 - - — the effect of Coulomb inter- 
action and the properties of the dot spectrum were taken 
into account. 

In the field of pumping, spin-dependent charge trans- 
port and spin transport have been studied extensively in 
the high-frequency rcgimei 27 ' 47 Adiabatic spin pumping 
has been explored in both the presence of spin-orbit cou- 
pling^ and finite magnetic fields* 4 ^ as well as for a quan- 
tum dot attached to ferromagnetic contacts^ Recently 



adiabatic spin pumping in a magnetic wire with domain 
walls was investigated^ Adiabatic pumping through 
double dots in the presence of spin-orbit coupling can be 
used to study the spin dynamics in a Pauli-blockade con- 
figuration^ Also charge and spin in ferromagnetic hy- 
brid structures pumped by magnetization dynamics have 
been investigated in detail^ An experimental realization 
of a spin pump used the modulation of a single quantum 
dot in the presence of a Zeeman fields 
In this paper we study a double-dot system with onsite 
Coulomb interaction contacted weakly to electronic reser- 
voirs. We consider the situations where both leads are 
normal conducting (N-DD-N), see Fig.[TJ where only one 
of the reservoirs is normal conducting while the other is 
replaced by a ferromagnetic contact (N-DD-F), and fi- 
nally the case where both contacts are ferromagnetic (F- 
DD-F). We investigate the effects of quantum charge fluc- 
tuations and the impact of different effective coupling to 
the hybridized double-dot states on the pumped charge. 
Based on these effects the pumped charge shows charac- 
teristic sign changes. 

Also in the presence of spin-polarized leads, these two ef- 
fects matter. With respect to an earlier work on pumping 
through a single interacting quantum dot^S the present 
study of a double dot reveals a number of strikingly dif- 
ferent effects, such as pure spin currents and the inversion 
of the transport direction due to polarized leads only. 
We use a real-time diagrammatic approach^ extended 
to the adiabatic regime^ to calculate the pumped charge 
and spin, taking into account Coulomb interactions. For 
the interdot hopping we treat both the case of weak and 
strong couplings. 

This paper is organized as follows: we introduce the 
model and the formalism used for our theoretical investi- 
gations in Sec. |lll In the following we present our results 
for charge pumping through a double dot coupled to nor- 
mal conductor leads (Sec. lIIl) ). We discuss spin pumping 
when one of the leads is spin polarized in Sec. IIV1 In 
the last part we consider the transport behavior in the 
presence of two ferromagnetic leads with arbitrary polar- 
ization angle (Sec. |VJ). We set h = c = 1 for the rest of 
this paper. 



II. MODEL AND FORMALISM 



Hamiltonian 



We consider two single-level spin-degenerate quantum 
dots coupled in series to each other and tunnel coupled 
to a left (L) and a right (R) lead, as depicted in Fig. [1] 
The Hamiltonian of this system is given by 
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FIG. 1. Schematic picture of the system. The left and right 
dot are coupled to the electron reservoirs L and R, with equal 
electrochemical potentials, by barriers with tunnel coupling 
strength Fl and Fr, respectively. The coupling between the 
single-level quantum dots is given by A and the two dots have 
time-dependent energy levels el if) and er (t). 



The Hamiltonian of the double dot is explicitly time de- 
pendent, 

H d d (t) = ^ Cq (*) n ° + Un L n R + u ' ( n Lf«i4 + "RT^iu) 

a=L,R 

-7$] ( d L d R° + 4.<A° 



a=t,i 



(2) 



The single level e a (t) of each dot a = L, R is modu- 
lated by periodically time-dependent gate voltages, with 
a frequency $7 and a respective phase shift. We de- 
fine the number operator for spin a =f, 4- 011 dot a as 
n aa = (v aa daa and the total number operator of dot a 
as n a = J2cr n acn where cr a(7 (d a<7 ) creates (annihilates) 
an electron with spin er on the dot a. Electrons on each 
of the dots are subject to on-site Coulomb interaction U' 
and to Coulomb interaction with electrons on the neigh- 
boring dot U. The charging energies are related to the ca- 
pacitances of the dots and are taken into account within 
a constant interaction model, see, for example, Ref. 154 . 
Hopping from one dot to the other occurs with the in- 
terdot coupling amplitude —A/2, where A is taken to be 
real and positive. It is useful to define the mean double- 
dot energy E = E (t) = [e L (£) + e R (*)] /2 =: E + 5E{t) 
as well as the level difference e = e (t) — el (t) — £r (t) =: 
e + 8e(t), separating each quantity into a time-averaged 
and a time-dependent part. We will later apply the 
convention of omitting the time argument and denoting 
time-averaged quantities by a bar also for functions con- 
taining these parameters. In the following we take U' to 
be the largest energy scale, and thus every dot is at most 
singly occupied, excluding exchange coupling. 

The non-interacting electrons in the left and right leads 
are described by 



H a =Yj tkci ka c ak a 



(3) 



H = H dd (t)+ J2 H * 



H< 



a=L,R 



mc i . (1) with the creation (annihilation) operators c akcr (c a k a ) of 

an electron with spin a and momentum k in lead a. We 



will include both cases of either ferromagnetic or nor- 
mal conductor leads. We assume the leads to have equal 
electrochemical potentials, p^ = (Ur, not influenced by 
the time-dependent gate voltages. The effect of a time- 
dependent bias is discussed in Ref . [22]. The tunnel Hamil- 
tonian describing the dot-lead coupling is given as 



-fftunncl — 2_^i 7 i ^« 
a— L,R h.a 



,t 



H.c 



(4) 



with the tunnel matrix elements Vl and Vr, which do 
not depend on momentum or spin. The tunnel rates are 
defined as r a = 2ir \V a \ p a for a = L, R. We define Y as 
the sum of the left and right tunnel rates, T = Tl + Tr. 
We assume all bare tunneling rates T, A to be time inde- 
pendent. We will see later on that the tunneling rates to 
hybrid dot states become effectively time-dependent via 
the time-dependent dot energy levels. Both reservoirs are 
taken in the wideband limit where their densities of states 
p a are energy independent. We will later consider ferro- 
magnetic leads by introducing spin-dependent densities 
of states in the respective lead a = L, R. The technical 
details for the treatment of ferromagnetic leads will be 
discussed in Sec. EI 



B. Real-time diagrammatic approach 

We describe the system by the reduced density ma- 
trix of the double dot; the degrees of freedom of the 
leads, taking up the role of baths, are traced out. We 
write the accessible states on the double dot as |x). The 
time evolution of the density matrix P* 2 = (x2|pdd|xi) 
is captured in the generalized master equation (kinetic 
equation) 



± P X2 

dt X1 



(t) 



>E 

Xi>X2 



rX2,X 2 
Xi.x'i 



(i)P^(t) 



y f dt'w x2 ^{tX)p x ht r ) 

Z^ / xi,Xi v ' J Xi v J 

, / . / J — oo 



(5) 



xi,x 2 



The first part represents the interdot dynamics via a 
Liouville superoperator, where we define the matrix el- 
ements (x2\L(t)pdd(t)\xi) = (X2|[#dd(*),Pdd(*)]|Xi) = 
Y^> ^ L X2 ' X ? (t)P x ? (t). The kernel W X2,X '? in the scc- 

^Xi.X 2 Xl.Xi w Xi v ' Xl.Xi 

ond part of Eq. ([5]) incorporates the lead-dot tunneling 
transitions between states x'i an d X2 a ^ time t' , and states 
Xi and Xi at time t and is calculated using a real-time 
diagrammatic approach^ We write the diagonal and off- 
diagonal elements of the density matrix in vector form, 

P = I {Pg } X = X ' , {P x } x #x' ) j enabling a matrix repre- 
sentation for the kernel W, and likewise for the Liouvil- 
lian L. The generalized master equation, Eq. ([5]), con- 
tains the full time dependence of the problem, which is 
affected by the time-dependent dot energy levels. In the 
limit when the lifetime of the double-dot states is much 



larger than the time scale given by the variation of the 
parameters, we can perform an adiabatic approximation 
along the lines of Ref. [23[ One first carries out a Tay- 
lor expansion, P (£') = £\ (*' - t) j /jl^jP (t), around 
time t of the matrix-elements of the reduced density ma- 
trix in the integrand of Eq. (|5|), accounting for a finite 
memory of the kernel. In addition, the internal time 
dependence of the kernel via the time-dependent param- 
eters is taken into account up to first order in $7. The 
adiabatic expansion of the kernel contains an instanta- 
neous part (i), where all parameters are fixed to their 
value at time t, X(t) —> X(t). It further contains a first- 
order correction term (a) which is found by considering 
X (r) ->■ X (t) + (t — t) ^r I T =t in the evaluation of the ker- 
nel and by taking into account systematically first-order 
terms in the time derivatives. This yields the two terms 
of the expansion 



w (t, t') -> wr (t - 1') + wr (t - f) 



(a) 



(6) 



The resulting parametric time dependence is indicated 
by the subscript t. With this, a set of equations for the 
instantaneous and the adiabatic part, P t z and P\. , of 
the reduced density matrix is constructed. The instanta- 
neous contribution of the density matrix fulfills the sta- 
tionary generalized master equation with system param- 
eters frozen to values at time t. The adiabatic contribu- 
tion accounts for the fact that the actual density matrix 
slightly lags behind its instantaneous value. 

We here restrict ourselves to the case of weak coupling. 
Therefore, on top of the adiabatic expansion for small Q 
we perform a perturbation expansion in the tunnel cou- 
pling between lead and dots. We take into account only 
terms up to first order in the tunnel coupling, which is 
valid if the broadening due to tunneling is smaller than 
the temperature broadening, T < k^T, where fee is the 
Boltzmann constant. We perform an expansion of the 
generalized master equation order by order in T, as used 
in Ref. [53| and subsequent works. This systematic ex- 
pansion ensures that all contributions to a given order 
are taken into account (it may fail in certain situations 
when higher order terms in Y play a role^ which is not 
the case in the regime considered here) . We will also con- 
sider the case where Y is not the only small parameter but 
the intra-dot coupling A and the level difference e have 
the same magnitude. Then coherences become important 
already in first order Y, and the perturbation expansion 
has to be performed taking into account consistently all 
three small parameters, see Sees. MI "SI and HV A[ as well 
as Appendix [K\ Gathering the terms of the lowest or- 
der in £1 and the small parameters of the perturbative 
expansion, we end up with 



0= W 



(».i) _ iL ( l ^)\ dM) 



(i,0) 



dt' 



W 



(*.l) ,;r(M)\ r>(a,-l) 



iL 



(7a) 
(7b) 



Here we have introduced the zero-frequency Laplace 
transform wf ;) = f*L 00 dt'W < P (t - 1'). The order in 



the perturbation expansion is indicated by numbers in 
the superscript of the respective quantities, i.e., first or- 
der in the transition elements W'*' , and the superscript 
of L^ 1,1 ' indicates that the double-dot Liouvillian enters 
this equation in the same order of the Kernel and that its 
time dependence is instantaneous. The instantaneous as 
well as the adiabatic reduced density-matrix elements are 
P (i ' 0) and p( a ~ l \ Note that the adiabatic correction to 
the reduced density matrix starts in minus first order in 
T; this is justified as these terms are proportional to fi/r, 
which is a small parameter in the adiabatic regime. Since 
we do not consider higher orders in the perturbation ex- 
pansion, the corresponding superscript for the above dis- 
cussed quantities is unambiguous, and is suppressed in 
the following for simplicity. 

Equations (|7a[) and (|7bj) have a very similar structure 
and therefore one can obtain Eq. (|7a|) from Eq. (|7b[) by 
setting the left-hand side to zero and replacing the re- 
duced density matrix by the instantaneous contribution. 
This will be used throughout the remainder of this paper. 

A similar expansion is performed for the current 
through the system. Generally, we can give the expres- 
sion for the pumping current through the double quan- 
tum dot as 



h(t) 



ee 



dt'W L (t,t') P [t') 



(8) 



with the unity charge e. The trace of all contribut- 
ing elements is performed by applying the vector e T = 
(1, . . . , 1, 0, . . . , 0), where the number of ones is equal to 
the dimensionality of the reduced Hilbert space. The 
kernel W (t,t') with the zero-frequency Laplace trans- 
form W takes into account processes which involve an 
exchange of charge between the left lead and the quan- 
tum dot subsystem, see Ref. [5j. In the case of zero bias, 
considered in this paper, the instantaneous contribution 
to the current, Jw, vanishes at all times and the only 
remaining term is the adiabatic contribution, 7™. In 
lowest (zeroth) order in the tunnel coupling we obtain 



I L (t):=l[ a) (t) 



ee W t P; 



(9) 



where we drop the superscript (a) for simplicity. Expres- 
sions §8§ and © represent the tunneling current into 
the left lead. Note that due to the time dependence, 
additionally also displacement currents are present. In 
Ref. [56| this has been studied in detail for a single dot. 
Importantly, the displacement current is a pure ac cur- 
rent and it averages out when integrating over one period 
of the pumping cycle T = 2tt/CI. In the following we are 
mainly interested in the pumped charge per cycle, hence 
only the tunneling current plays a role and we discard 
the displacement current from now on. The number of 
pumped electrons is then given by 



1 

N = - 



T 



We consider the system in a time-dependent steady state, 
i.e., no transient behavior; therefore the total charge of 
the double dot is conserved after one pumping cycle.— 



III. NORMAL CONDUCTOR RESERVOIRS 

In this section we discuss the pumped charge through 
the double-dot system in contact with two normal- 
conducting leads (N-DD-N). 

A. Weak interdot coupling 

We first consider the situation where the two dots are 
weakly coupled to each other, A < k-gT, and the energy 
difference between a singly occupied left and a singly oc- 
cupied right dot is small as well, e < kftT : such that both 
parameters are on the order of the coupling strength, 
F ~ A ~ e. We account for expressions in the mas- 
ter equation up to first order in these parameters. In 
this case a rigorous expansion in the coupling param- 
eters (up to lowest order in both lead-dot and dot-dot 
couplings) and the level difference e is performed. This 
system has been considered in the static case in Ref. l29i . 
For the remainder of this section we assume that not 
only U' but also the charging energy U is much larger 
than all other energy scales (such as temperature, the 
modulation frequency, and the level difference e). There- 
fore the double-dot system can only be singly occupied 
or empty. In this case the states |L<r) and |Rcr), with the 
electron with spin a =f, 4- m the L or R dot, are almost- 
degenerate quasi eigenstatcs of the system and coherent 
superpositions of |Lcr) and |Rc) play an important role, 
see Appendix El These coherent superpositions are cap- 
tured in the off-diagonal elements of the reduced density 
matrix of the double-dot system. Even in the absence of 
bias, they do not vanish for an asymmetrically coupled 
double dot due to the time dependence of the system. 
The vector of the reduced density-matrix elements is 
P = (P ,P L ^P Li7 P Rt7 P Ri ,P^,P^,P^lP^), where 
we write the diagonal elements of the density matrix 
P* =: P x . We end up with a master equation for the 
occupation probabilities for an empty dot Pq and for 
a singly occupied dot Pi = Pl + Pr, where the total 
occupation with different spins is Pl = P^ + Plj, an d 
Pn = PR t + PRl- We find 



dt 



P = -217+ (E) P + c z XT.r (E) S 

+ir/-(p)p, 



111) 



dth (t) 



(10) 



and Pi is obtained via the probability conservation P W + 

Pf = 1 and P + -P{ = 0, for the instantaneous 
and the adiabatic parts of the reduced density-matrix 
elements. The tunnel coupling asymmetry is given by 
A = (Tl — Tr) /r taking values between —1 and 1, where 



zero is the case of symmetric coupling to the leads. The 
Fermi function is /+ (E) = efl i +1 and /~ (E) = 1 — 
/ + (E) with the inverse temperature j3 = l/k&T. The 
vector e z projects out the z component of the pseudospin 
vector which captures the off-diagonal elements and the 
difference in the occupation of the left and the right dots, 



(12) 



The dynamics of the pseudospin is described by a Bloch- 
type equation 



j/ = c z Ar (f+ (E) P Q - if- (E) P 1 




1 



(13) 



Tf-(E)S + BxS 



The time evolution of the pseudospin has a contribution 
due to the accumulation of the pseudospin z component, 
i.e., an unbalance of the left-right occupation, cf. first 
line of Eq. (|13j) . Accumulation occurs only for asymmet- 
ric lead coupling. Relaxation of the pseudospin, given by 
the first contribution of the second line of Eq. ([13]) takes 
place independently of the coupling asymmetry. Further- 
more, we find a precession of the pseudospin around an 
effective magnetic field, via which the internal dynamics 
of the double dot enter. It is given by 



-A 



B 



(14) 



and intra-dot Coulomb interaction is found in Ref. l29|. 
The effective field due to the interdot tunneling lies in 
the x, y-plane which is spanned by the pseudospin con- 
tribution due to the coherent superpositions of |Lcr) and 
|Ra). 

The instantaneous solution of the master equation, 
Eqs. (jlip and (fi"3)) , is given by Boltzmann distribu- 
tions. The average instantaneous occupation number on 
the double dot (n) W = • Pq + 1 • P\ is given by 



,(<> 



l + 4e -/3B 



(17) 



The factor four stems from the total fourfold degeneracy 
of the singly occupied state. For the same reason, we also 
have 



S {%) (t) = 



(18) 



in lowest order in the tunnel couplings T and A and the 
level difference e (this holds as long as no bias voltage 
is present). As opposed to the instantaneous solution, 
the adiabatic correction of the pseudospin S^' -1 ) does 
in general not vanish. This is due to an occupation dif- 
ference in pseudospin space, introduced by the time de- 
pendence and asymmetric coupling to the left and right 
leads, A ^ 0. Therefore, in the case of symmetric dot-lead 
coupling, the adiabatic correction for the pseudospin also 
vanishes. For the adiabatic current we find from Eq. ([§]), 



ed (n) 
Ih{t) -2—t 



d) 



(19) 



We note that instead of the bare level difference e, a 
renormalized one is entering the effective field, given by 



AA 2 



e + X Y J(E) 



(15) 



This expression is explicitly time dependent via e and E. 
The function <fr entering the renormalized level difference 
is given by 



i (E) = Re 



1 2 2tt 



,i + <W), (16) 



where ip denotes the digamma function. The renormal- 
ization arises from quantum charge fluctuations. It enters 
the z component of the effective field, which acts as a Zee- 
man field for the pseudospin. This Zeeman field affects 
the pseudospin dynamics through a precession around 
the z axis, thereby coupling to the x and y component 
of the pseudospin, which arise due to coherent super- 
positions of | Lit) and |R<r). The level renormalization 
is a pure Coulomb interaction effect, i.e., it vanishes for 
U = U' = 0. Furthermore, the renormalization is zero for 
symmetric coupling A = 0. Here, Coulomb interactions 
are large and e cu toff provides the cutoff energy. The full 
expression for the renormalization of e for finite inter- 



(i - A 2 ) (ir 2 [/- (£)r 



Note that the parameters e and E depend on time and 
that the pumping current is proportional to the time 
derivative of the occupation number, depending on E. 
Therefore, a necessary prerequisite for a nonzero adia- 
batic current is an explicitly time-dependent E. Impor- 
tantly, also the prefactor of -4 (n) is in general time 
dependent. This is the necessary condition for a non- 
vanishing average pumped charge.— 
If A = we find that the time-dependent current through 
each contact a = L,R is given by I a (t) = §Jj (n) . 
Therefore the currents injected into the left and the right 
leads are equal at any time and the pumped charge per 
cycle vanishes; it is therefore directly sensitive to the cou- 
pling asymmetry. 

The pumped charge is obtained from Eq. (|TU|) . We are 
interested in the regime of bilinear response for the mod- 
ulation, which is valid if the pumped charge per infinites- 
imal area in parameter space varies little within the area 
enclosed by the cycle. Then the pumped charge is pro- 
portional to the cycle area, A = L dtSeSE. We find 
the number of pumped charges N per area in parameter 




1.5 r 



FIG. 2. (Color online.) Plot of the pumped charge JV with 
respect to the mean bare level difference e for E = 1.5/cbT 
and asymmetric tunnel coupling, A = 1/2. The node N — 
occurs when e rcn = 0. The other parameters are A = T, 
k B T — 2F, e cuto ff = 50fc B T. 



space as a function of the average quantities e and E, 



N 

T 



-A(l-A 2 



A 2 



d(n 



d) 



(l-A 2 )([E/-(S)] 2 + 4 n )+A 2 



dE 



(20) 

The pumped charge has a peak when E is close to reso- 
nance and fulfills the relations 



N (e ro „) =-N (-e ren ) 
N (A) = - TV (-A) . 



(21a) 
(21b) 



Inverting both parameters consequently maps the func- 
tion N onto itself again. As discussed above, the pumped 
charge vanishes if A — 0, as well as at A = ±1 when one 
of the leads is completely decoupled. Since the pumped 
charge also vanishes if e ron = (see Fig. [2]), the level 
renormalization can be directly read out by means of 
pumping through the double dot, when scanning through 
the time-averaged left-right level difference e. This im- 
portant property occurs due to the following reason: the 
prefactor of 4? {n) l ' is an even function of the renormal- 
ized difference of the left and right level position e len (t); 
therefore in the limit of small pumping amplitudes (bilin- 
ear response, i.e., only the time dependence of e matters 
in the prefactor) a sign change in the average e ron has the 
same effect as a shift of the modulation by it. Doing the 
average over one pumping period, the transport direction 
is therefore reversed. 

The level renormalization due to quantum charge fluc- 
tuations, which can be measured in the pumped charge, 
is distinguishable from possible level renormalization ef- 
fects due to an energy-dependent density of states by its 
temperature dependence, which is logarithmic for large 
k B T. 



B. Strong interdot coupling 

A finite interdot coupling amplitude A leads to a hy- 
bridization of the energy levels of the right and the left 



dots. In the limit of strong coupling between the two 
dots, A > r, we write the single-particle part of the 
double-dot Hamiltonian in terms of bonding |b) and an- 
tibonding |a) states, '^2, r . a e v (£) d^d^, where r\ = b, a 
and a =f, \.. The energy of the bonding and antibonding 
state is given by 



£b/a 



et\^ 2 



(22) 



The corresponding operators creating (annihilating) an 
electron in the bonding and the antibonding states are 
related to the ones creating an electron in the left or the 
right dot by the equality 



U L/R 
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V2 



1± 



Va 2 
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(23) 

where we omitted the spin indices for simplicity. The 
bonding |b) and antibonding |a) states have an energy 
difference much larger than the level broadening. This 
means that the interdot dynamics is much faster than the 
dot-lead hopping. Therefore, coherent superpositions of 
bonding and antibonding states are suppressed in low- 
est order in the tunneling and no dynamics of the off- 
diagonal elements of the reduced density matrix have to 
be considered (in contrast to the previous case discussed 
in Sec. IIII A[) . We now take U to be finite and we there- 
fore also consider the occupation of the left and right 
dots with one electron each. We still assume large on- 
site Coulomb interaction, inhibiting double occupation of 
each single dot. As a consequence, the doubly occupied 
(fourfold degenerate) states are |LerRcr'). We can now 
write the double-dot Hamiltonian ^ E x \\) (x| in the ba- 
sis of the eigenstates |x)- Here the eigenenergy E of the 
empty dot |0) equals zero, the eigenenergy £"b of the spin- 
degenerate bonding states |ber) = c^JO) is e^ and anal- 
ogously the eigenenergy Eg, of the spin-degenerate anti- 
bonding states \&a) = (4 a |0) is e a . Finally the four dou- 
bly occupied states \LcrRa'), with a =t,| and a' =t,4, 
have the eigenenergies E aa i = 2E + U. The tunnel cou- 
pling between the double dot and the leads a = L, R is 
captured via effective rates for tunneling through the hy- 
brid states |b) and |a). These rates are explicitly time 
dependent and are given by 



arj 



VA 2 



(24) 



Tunneling through the hybrid single-particle states is de- 
noted by the subscript r\ = b, a. The difference of the cou- 
pling strengths for the two transport channels is propor- 
tional to a factor of e/VA 2 + e 2 . (This factor is related 
to the relative position of the bonding and antibonding 
states with respect to the localized states L and R.) In 
order to bring out this important property we use the 
following notation: if i] is used as a variable rather than 
a coefficient it takes the value +1 (—1) for b (a); equally 
if a = L, R is used as a variable it takes the values +1 for 
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FIG. 3. Plot of the pumped charge N as a function of 
E for different values for the Coulomb interaction U — 
0, 5/cbT, 10/cbT (dashed, dotted, solid), in the symmetric case 
e = 0, A = 0, and the interdot coupling is A = k B T and 
k B T = 2I\ 



L and —1 for R. The sum of these rates is T v = ^ Q T ari . 
With this, the master equation for the occupation prob- 
abilities P , P n = P, n + P r)i , and P d = J2 a ,a'=t,i -Fto-Rff' 
reads 

yPo= Y, (-2T n f + (e ri )P + T n f-(e n )P n ) , (25a) 



dt 



dt 



77= b,a 



P„= 2rj + (e v )P + (-TJ-(e 11 ) 



(25b) 



- 2I>/+ (e rl + U) 1 P n + T n f- (e n + U) P d , 

(i) 

and P c \ is obtained via the probability conservation P + 

pW + pW + pW = ! and p^«) + pM + pW + p< (a) = Q 

The instantaneous probabilities arc again given by the 
Boltzmann distribution. 

As described before, in Sec. Ill Bt the instantaneous oc- 
cupation probabilities and their first-order corrections in 
the adiabatic expansion are found as solutions of this 
master equation. The adiabatic current is then calcu- 
lated straightforwardly as 



M'l-E^f-^' + ip, 



7]— b,a 



(26) 



The current expression is divided into two different con- 
tributions, Ih (t) = Ih,b (t) + Ih,a (*)• It is apparent that 
the dot state transitions |0) <-> |b) and |a) <-> |d) occur 
with the rate Tb. And likewise, the transitions |0) <-> |a) 
and |b) O |d) are due to tunneling with the rate T a . 

In bilinear response for small pumping amplitudes, the 
pumped charge through the double dot becomes 



A ^ dl { r„ ) dE 



p(i) 
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p. 
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(27) 



The derivative with respect to E of the sum of the prob- 
abilities is always negative for both the bonding and the 
antibonding contributions. In contrast, the prefactor, 
namely, the derivative with respect to e of the relative 
effective coupling, has opposite signs for the different hy- 
brid states, since the time-dependent part of the effective 



coupling comes always with opposite signs, respectively, 
see Eq. (|23|) . 

As a direct result, one finds that the terms N\, and 
A a , derived from the current separated into iL,b (t) and 
Il,& (t) , always contribute to the pumped charge per area 
with opposite signs. This result is a central point of our 
paper; this double-dot feature cannot be provided by a 
single-dot system and is of fundamental importance for 
many of the effects discussed in the following. 

In Fig. [3] we plot the pumped charge as a function of 
the mean level E in the case of spatial symmetry e = 0, 
A = and for different values of U. Importantly, the 
finite Coulomb interaction shifts the resonance positions 
and furthermore enhances the pumped charge. The two 
resonant peaks have opposite signs; they appear when the 
addition energy to go from an empty dot to the ground 
state of the singly occupied dot is at resonance, £b ~ 
and when the addition energy to go from the singly oc- 
cupied ground state to a doubly occupied dot is at reso- 
nance, f a + [/«0, except for the common temperature- 
dependent shift due to different charging and decharging 
rates. Since the main contribution at eb ~ comes from 
tunneling involving the hybrid state |b) and the main 
contribution at e a + U ra comes from tunneling involv- 
ing the hybrid state |a), the pumped charge has opposite 
signs at the two resonances. The pumped charge from the 
two contributions does not cancel out for U — because 
of the finite level spacing. This result agrees with Refs.[2J 
andl59l. where charge pumping through a double dot with 
nonintcracting spinless electrons was considered. 

If the double-dot parameters arc such that the level 
splitting |eb — £a| is much larger than k^T, the probability 
of having state |a) occupied vanishes. Thus, the trans- 
port close to the single-electron resonance occurs solely 
via charging and decharging the single-electron ground 
state |b), via the processes |0) -O- |b), contributing with 
rate Tb- Equally, at the resonance for the transition be- 
tween singly and doubly occupied double dots, only the 
transitions |b) •f* |d) contribute with rate F a . 

In the regime when the excited state is still thermally 
accessible, |eb — e a | ^ k^T, also the transitions with the 
excited state, |0) «-» |a) and |a) <H> |d), start to con- 
tribute. By introducing a spatial asymmetry to the dou- 
ble dot, choosing e / or A ^ 0, the two contribu- 
tions (for transport involving the ground state and trans- 
port involving the thermally excited state) are changed 
by different amounts. In spite of the fact that the 
larger time-dependent current flow stems from transport 
through the ground state, for the pumped charge, i.e., 
the time-avcraged quantity, the contribution through the 
ground state can become the minor one. This means that 
by choosing appropriate double-dot parameters charging 
and decharging of the ground state leads to smaller di- 
rected transport than charging and decharging of the ex- 
cited state, a unique feature of time-averaged transport 
due to time-dependent fields. Thereby, we can achieve 
a sign change in the pumped charge at one of the res- 
onances, see Fig. SJa), (analogously for the other reso- 
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FIG. 4. (Color online.) (a) Plot of TV as a function of E for 
the asymmetric case (blue dashed) and the symmetric case 
taken from Fig. [3] as reference (black solid), (b) Map of the 
peak height of the pumped charge N at the single-electron 
resonance (dashed vertical line in (a)) in dependence of e and 
A = (T-l — Tr) /r. The empty and filled circle indicate the 
symmetry configurations taken in (a). For all plots the in- 
terdot coupling is A = k B T and the Coulomb interaction is 
U = 10fc B T. The thermal energy is k B T = 2T. 



nance at e a + U ~ a sign change appears for A — > — A 
or e — > — e). 

We plot the peak height of the pumped charge close to 
the resonance l^ ~ for the full parameter spectrum of A 
and e in Fig.UJb). There, we see that for most parameter 
configurations, the bonding channel is dominant hence 
the peak is positive. The blue regions (brighter areas in 
the lower right and the upper left) show the parameter 
regimes where a sign change appears due to the dominant 
transport through the thermally accessible antibonding 
state. These sign changes in the peaks with changing 
asymmetry will however vanish, once the level splitting 
becomes significantly larger than k^T and the excited 
state is no longer thermally accessible. We also find that 
the pumped charge is point symmetric with respect to 
the parameter A and the bare, time-averaged e, N (A, e) = 
N(-X,-e). 



IV. SPIN PUMPING THROUGH N-DD-F 

We now replace one of the contacts by a ferromag- 
netic lead (we choose the right one), thereby breaking 
spin-rotation invariance which enables spin in addition 
to charge pumping. Spin and charge pumping through a 
single interacting quantum dot in the presence of a fer- 



romagnetic lead has been studied before^** spin pumping 
through a double dot is particularly promising due to the 
following reason: the possible sign reversal of the charge 
transport which in the weak-coupling regime is taking 
place at the renormalized level difference e rcn being zero, 
and which in the strong-coupling case is due to thermal 
accessibility of the excited level, is expected to affect the 
transport for different spin differently 
We now discuss the representative situation at the res- 
onance between empty and singly occupied double dots 
and therefore restrict the calculation to infinite Coulomb 
interactions U and U' for the remainder of this paper. 
In a spin-polarized contact a the density of states is spin 
dependent, p a ^ ^ p Q j_. The spin polarization strength of 
lead a is defined as 



Pa-\ - Pal 
Pat + Pal 



(28) 



Due to the spin-dependent density of states also the tun- 
neling rates to the right lead become spin dependent 



(1 + <rp a )T 



(29) 



where a = +1 for r a -t- and a 



The total tunneling rate for lead a is T a = ^ (Tai + r Q j_). 



-1 for T a i and a = L, R. 

r 

In this section only pr is different from zero hence 
Tlct = Tl- We choose the polarization axis of the spin 
in all parts of the N-DD-F system along the axis of the 
majority spin of the ferromagnetic contact and we can 
study the dynamics of electrons with spin up and spin 
down separately. Due to the spin-dependent tunneling 
rates, the spin-up and the spin-down channels have dif- 
ferent pumping dynamics, leading to a generally nonzero 
net spin transport. Analogously to Sec. IIH( wc discuss 
in the following the two limits for weak and strong in- 
tcrdot coupling. We find that the equations for the spin- 
resolved pumped charge are formally equivalent to the 
results of the unpolarized case, see Sec. IIII1 where one 
replaces T a — > T aa . We will show that the two pre- 
viously discussed different regimes both lead to a pure 
pumped spin current, relying on different effects. 



A. Weak interdot coupling 

We first discuss the case of weak interdot coupling, 
where A ~ T ~ e, analogous to Sec. IIII Al In order to 
calculate the total pumped charge N = Nf + JVj, and 
the total pumped spin N s = Nf — TVj,, we evaluate the 
number of pumped electrons with spin a. We obtain an 
expression similar to the unpolarized case from Eq. (|20l) , 
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FIG. 5. (Color online.) Plot of the pumped charge through 
the up spin state Nf and down spin state iVj, (upper plot), 
and the total pumped charge A*' = Nf + N± as well as the 
total pumped spin N — Nf — N± (lower plot), for different 
coupling symmetries (a) A = —0.6 and (b) A = —0.1. The 
polarization strength of the right reservoir is pr = 0.5. The 
other parameters are A = T, E = 1.5kBT and feT = 2I\ 



with the spin-dependent quantities 



r 



Iia 



e re 



Tlct + Tro- 

€+ 2, 



4>(E) 



(31a) 
(31b) 



Importantly, the coupling asymmetry is spin dependent 
here, more specifically, for the spin-dependent coupling 
asymmetries we always find A-j- < A < A^. This also 
makes the rcnormalized energy-level distance e ren ,<T spin 
dependent. As we have elaborated in Sec. IIII Al there 
are nodes in the pumped charge whenever the left and 
right tunnel rates are equal, and when the time-averaged 
renormalizcd level difference e rcn is zero. In analogy, we 
find nodes for the spin- resolved pumped charge N a : zero 
net transfer of particles with spin a occurs if A CT = and 

if Cren^cr = 0. 

The spin-resolved pumped charges Nf and N± are plot- 
ted in the upper panels of Fig. [5] for different values of 
the coupling asymmetry. Clearly, the number of pumped 
charges with spin up and down differs strongly from each 
other leading to a finite spin transport. If the coupling 
asymmetry is such that Tr^ < Tl < Tr-)-, then spin- 
up and spin-down electrons are even mostly pumped in 
opposite directions as is shown in the upper panel of 
Fig.E(b). 

The lower panels of Fig. [5] show the total pumped 
charge and the total pumped spin. We see that the 
pumped charge and spin have different nodes: whenever 
Nf = ATj, the pumped spin vanishes, as indicated by the 
vertical dotted line in Fig. [3] (a), while the pumped charge 
can still be finite. On the other hand, when Nf = —N± 
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FIG. 6. (Color online.) (a) Plot of the pumped charge re- 
solved for all four possible channels N vcr , r\ = b, a and a =|, F, 
(b) Spin-resolved pumped charge, Nf and N±. (c) The total 
pumped charge N — Nf + Ni as well as the total pumped spin 
N = Nf— JVj.. For all plots the coupling symmetry is A = 0.6, 
the polarization strength of the right reservoir is pn = 0.5, 
and the other parameters are A = faT, E — 1.5fcBT and 
k B T = 2r. 



the pumped charge vanishes, as indicated by the vertical 
dotted line in Fig.[5](b), while the pumped spin can still 
be finite. This results in both the possibility to pump 
charge in the absence of net spin transport, and quite 
more intriguing, the possibility to pump spin without net 
charge transport by electrical control only. In the regime 
studied here, pure spin pumping in the absence of charge 
pumping is due to spin-dependent quantum charge fluc- 
tuations, induced by Coulomb interactions, resulting in 
differing spin-dependent level rcnormalization effects. 



B. Strong inter dot coupling 

We now turn to the strong interdot coupling regime, 
A > T, which is studied in Sec. IIII Bl for unpolarized 
leads. Similar to Eq. (|27[) . we obtain for the spin-resolved 
pumped charge through the two states r\ = b, a, 



N, 



IjCT 



A 



ld_ 
2dl 



Ln 



0PrT 



Ri/ 



dP, 



(<) 



dE 



(32) 



These quantities depend explicitly on the spin and the 
state r\ through which the charge is pumped and therefore 
all four quantities are in general expected to be different, 
leading to possible pure dc charge as well as pure dc spin 
transport, sec Fig. Htc). 

To explain the origin of different Nf and iVj, we con- 
sider Figs. [SJa) and HJb). In Fig. EJa) we resolve the 
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FIG. 7. A schematic picture of the Hamiltonian with ferro- 
magnetic leads. The magnetization of the leads is given by 
p L and p R . 



pumped charge for each channel separately, namely, for 
a spin-up electron through both hybrid states, iVbf and 
A^^, and likewise for a spin-down electron, Afbi and Ng±. 
For the plots in Fig. [5] we take et> ~ in the case of 
strong coupling asymmetry A = 0.6 where we consider A 
and e such that the excited state (antibonding) is ther- 
mally accessible. Figure [6ja) shows contributions due to 
transport through the bonding state (red), which have 
always an opposite sign to the charge pumped through 
the antibonding state (blue) . Since the spin polarization 
enters Eq. ([24]) only as a modification of the prefactor Tr , 
the important sign dependence of the time-dependent 
part of the effective tunnel coupling for different hybrid 
channels is not altered. Consequently as discussed in 
Sec. IIII B[ also the spin-resolved pumped charge has op- 
posite signs for different hybrid channels. Whether the 
charge pumped through the excited, antibonding level 
can dominate over the contribution through the bond- 
ing state, depends strongly on the coupling asymmetry 
to the leads. In the example shown here, for the spin- 
down channel the increased effective coupling asymmetry 
allows for dominant transport through the antibonding 
hybrid state. This leads to a sign change in the pumped 
charge with spin down N± as a function of e, in the limit 
of a thermally accessible excited state, see Fig. EJb) . In 
contrast, there is no configuration in which the transport 
of electrons with spin up takes place prevalently through 
the excited state since A-j- < A is too small to invert the 
transport direction. The important result is that N-^ and 
A^ can have opposite signs for certain values of e. There- 
fore we find again the possibility of pure spin pumping 
while the pumped charge vanishes and vice versa. 



NONCOLLINEAR FERROMAGNETIC 
RESERVOIRS 



In this section we discuss the pumping characteristics 
through the double-dot device coupled to two differently 
polarized ferromagnetic leads, F-DD-F, as sketched in 
Fig. [7J We concentrate on the strong intcrdot coupling 
limit A > r and start from the Hamiltonian discussed 



in Sec. Ill At it is now important to keep track of the 
different tunnel coupling for minority and majority spin 
of the different reservoirs. 

The polarization strength p ai Eq. (|28p . and the po- 
larization direction defining the polarization axes n a are 
different for the two leads a = L, R. We choose the spin 
quantization axis of each lead along the respective mag- 
netization axis, nL and /Tr, where the two vectors enclose 
an angle tp. Following Ref. |42|, we now take a coordinate 
system such that its basis vectors e x , e y , and e z align 
with n L + «r, «l — «r, and rtr, x tir, respectively. For the 
dots' spin quantization axis we choose the z axis of this 
coordinate system, which is perpendicular to the plane 
spanned by the leads' polarization vectors, see Fig. [7] In 
this basis, the left and right lead magnetization vectors, 
PL = Pl«l and p R = prtir , are given as 

/ cos ((p/2) \ / cos (<p/2) 

Pl=Pl sin (tp/2) , pr = pr - sin (ip/2) 

\ ° / V ° 

(33) 
It is useful to express the tunneling Hamiltonian, Eq. Q , 
in terms of creation (annihilation) operators for an elec- 
tron with a spin along the spin quantization axis of the 
respective system part. We therefore write 



H 



tunnel 



£ |E{i + ( e '" /4 ^t + e--% 



a=L : R 



+4_ (-e"^ 4 d Qt + e~ ia ^d ai ) } + H.c. , 

(34) 

where c ak _ i _ (c a k±) creates (annihilates) an electron with 
momentum k in lead a with majority/minority spin of 
lead a. Again, we use the notation that whenever a is 
used as a variable rather than a coefficient, it takes the 
values +1 for L and —1 for R. In this representation of 
the tunneling Hamiltonian it becomes apparent that the 
diagonal elements of the reduced density matrix couple 
to the off-diagonal ones with respect to the real spin. 
The reason for this coupling of diagonal and off-diagonal 
elements is that tracing out the ferromagnetic leads, we 
obtain dynamics for the reduced density matrix of the 
double dot which do not conserve spin. Therefore, off- 
diagonal elements of the density matrix have to be taken 
into account here, and we consider the density-matrix el- 

ements P = (p , P bt , P H , P at> P ai , P#, P$, P$, P^) . 

The off-diagonal elements and the difference of the oc- 
cupation number for different spins are contained in the 
vector 
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(35) 



'/! 



for 7/ = b, a. The generalized master equation for the 
elements of the reduced density matrix reads 
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Again p> is obtained via the probability conservation due to a spin dependence of interfacial phase shifts^ as 
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last line in Eq. (|36b[) represents a Bloch-type equation 
for the spin expectation value, which is affected by relax- 
ation and accumulation of spin and by a rotation due to 
an effective magnetic field which reads 



B 



q 7 Pa*- ar 

IT ' * 



(37) 



a=L,R 



The function </> is defined by Eq. (jTHJ) . This effective mag- 
netic field is induced by coherent transitions between spin 
states and is different for the bonding and the antibond- 
ing channels. Due to the absence of a bias voltage, we 
find that the instantaneous density-matrix elements do 
not differ from the previous case without magnetization. 
In particular, we find that the instantaneous spin expec- 
tation value is zero, Stf = 0. The adiabatic correction 
to the spin is different from zero; nevertheless, in the 
absence of bias, the adiabatic spin vector S,, ~ turns 
out to be parallel to the effective magnetic field^ hence 
no precession takes place and the effective magnetic field 
does not influence the transport dynamics. 

We calculate the pumping current through the double- 
dot system in the presence of arbitrarily polarized ferro- 
magnetic leads and find 
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with the definition of 
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(38) 



(39) 



The sign of the correction term to the pumping current 
with respect to the nonmagnetic case [second term in 
Eq. ([38|l ] differs depending on the difference in the lead 
polarizations. This means that the pumping current can 
be reduced in presence of differently polarized ferromag- 
netic leads, called the spin- valve effect ; 38 ' 39 or enhanced, 
showing an inverted spin-valve effect. Such an observa- 
tion was made for the pumping current through a single 
quantum dot attached to ferromagnetic leads i^ 3 - In this 
case of a single dot, the spin-valve effect for the pumped 
charge can only be inverted, when inducing a strong spa- 
tial asymmetry regarding the tunnel coupling. Consider- 
ing a static transport bias, an inverted spin-valve effect 
has, for example, been found for a single quantum dot 




well as for a double quantum dot due to charge fluctua- 
tions in the nonlinear bias regime4£ 
Qualitatively very different features for the spin- valve ef- 
fect are found when pumping through a double dot. We 
show in the following that this is due to effective spin- 
dependent asymmetries in the coupling to the different 
hybrid states. 

For the pumped charge per area of the pumping cycle (in 
bilinear response) we find 



(40) 



We use the abbreviation 7V P to indicate the pumped 
charge in presence of spin polarized leads to be con- 
trasted with the pumped charge N° = N(tp — 0) = 
N(pl = pn = 0), in the presence of normal conduct- 
ing leads [see Eq. (|27|) ]. According to Sec. MI 51 we sep- 
arate the two contributions corresponding to pumping 
through the bonding and antibonding states, contribut- 
ing with opposite signs to the pumped charge, N^/A > 
and N p /A < 0. The total pumped charge in the presence 
of polarized leads with respect to the pumped charge for 
vanishing polarization, N p /N°, depends on the strength 
of the polarization of the two leads and on their respec- 
tive angle tp. 

The quantity N p /N° is independent of E; it is shown 
for different parameter regimes in Fig. [S] As a function 
of the angle ip, the relative pumped charge N p /N° is 2ir 
periodic, and it equals one if tp = 0. The stronger the 
polarization of the leads, the more the functional depen- 
dence on tp deviates from a cosinelike behavior^ 9 - 
Three different parameter regimes are distinguished in 
panels (a)-(c), where we choose p\, = pn = p. The up- 
per plot, Fig. |U[a) shows results for a symmetric double 
dot, A = and e = 0, where the usual spin-valve ef- 
fect is found. The spin- valve effect is most dominant 
when the left and right polarization vectors are in an- 
tiparallcl alignment, tp = it. When the double dot is 
made asymmetric, the spin-valve effect can be inverted, 
N p > N°, ergo we find an anomalous spin- valve effect, 
see Fig. Hfb). In this case, the deviation from the co- 
sine behavior goes so far that for large polarizations, the 
maximum effect can even be found at angles away from 
7r. This implies that in the limit p — $■ 1, the relative 
pumped charge equals one for tp = and zero elsewhere. 
The most surprising result is that even the transport di- 
rection can be inverted as compared to the normal case, 
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FIG. 8. Plot of the ratio of the pumped charge for polar- 
ized leads to the pumped charge for non-magnetic contacts 
N p /N°, as a function of the magnetization angle tp for in- 
creasing polarization strength p — pl = Pr = {0.2,0.4,0.7}. 
(a) Plot for the symmetric case e = 0, A = 0. (b) Plot for 
the case of e = —ksT, A = 0. (c) Plot for the strongly asym- 
metric case e = — ksT, A = 0.6. For all three plots we take 
A = k B T, k B T = 2I\ 



leading to N p /N° < 0, sec Fig. [S^c). Also here, there 
is a shift of the maximum effect away from <p — -k as in 
Fig.[HIb), which becomes apparent for large polarizations 
and which is not shown in this figure. 



A. TMR of the pumped charge 

To classify these three cases in a more convenient way 
we define in analogy to the static cas o 38 i 46 a tunneling 
magnctoresistancc (TMR) for pumping 



TMR = 



N° -NP 
N° + NP' 



(41) 



Values between and 1 represent the case where trans- 
port is suppressed due to the spin-polarized leads, see 
the example in Fig. (Ha). Values between —1 and 
on the other hand indicate enhanced transport (reversed 
spin- valve effect), see the example in Fig. |S(b). In the 



more exotic case, where transport is reversed, the ab- 
solute value of the TMR is bigger than one and it di- 
verges when N p = -N°. The TMR is plotted as a 
function of the double-dot parameters in Fig. [51 where 
we show the case of antiparallel alignment of the mag- 
netizations of the two leads, tp = ir and p = p^ = Pk- 
The three different regimes are represented by a color 
scale: the red area denotes the ordinary spin- valve effect 
(0 < TMR < 1), and in the blue area the spin- valve ef- 
fect is reversed (— 1 < TMR < 0). The parameter area in 
which the charge transport itself is inverted (|TMR| > 1) 
is drawn in yellow, and the black line dividing the yellow 
area is indicating N p = —N°. For comparison, the pa- 
rameter choices for the three plots of Fig. [8] are indicated 
by the corresponding letters. 

The inverted spin- valve effect, — 1 < TMR < 0, can 
be easily explained for this case of (p = %. Now Nf, the 
number of pumped electrons with spin up, relates, e.g., 
to the majority spin in the left lead and the minority spin 
in the right lead. The shift between iV-j- and N±, similar 
to what is shown in Fig. [BJ depends on the polarization 
strength, in particular, the shift is zero for p = 0. This 
can lead to an increase or a decrease in the total pumped 
charge as a function of e and A, implying that the total 
pumped charge in the presence of polarized leads, N p , 
can be larger than for normal conducting leads, N°. This 
effect can occur independently of the accessibility of the 
excited (antibonding) state, as it does not rely on a sign 
change in iV-j- or N±. 

When |ejb — e a | i$ k^T parameter regimes exist where 
the transport changes direction, i.e., iV p and N° have dif- 
ferent signs [as depicted in Fig. Et c )] ■ As we have stated 
before, the polarization changes the transport through 
the spin-dependent bonding and antibonding channels 
differently. Consequently the transport through the ex- 
cited antibonding state becomes dominant at different 
parameter configurations, depending on the polarization 
of the leads. Therefore we can obtain different nodes 
N2 = (contours separating yellow and blue 
N p = (contours separat- 
ing yellow and red areas in Fig. \§§. In between these 
nodes, the charge transport is reversed. The regions of 
inverted transport can be enlarged by increasing polar- 
ization strengths p. 

We compare this with the tunneling magnctorcsistance 
for a double-dot system subject to a static bias in the 
linear-response regime. All parameters are therefore time 
independent. We define the static tunneling magnetore- 
sistance as 



for iV° 

areas in Fig. [9]) and N£ 



TMR, 



G° -G p 
G° + G p 



(42) 



via the linear conductance in the presence of normal 
leads, G°, and in the presence of spin-polarized leads, 
G p . Their explicit form is given in Appendix [Bj One 
can show that for all possible parameters in the regime 
of strong intcrdot coupling 



< TMR, 



<1 



(43) 
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FIG. 9. (Color online.) Overview map of the tunneling mag- 
netoresistance of the pumping system for p = 0.7, tp = n and 
A = ksT; the TMR is independent of E. Red represents 
the spin- valve effect, blue its inversion, and the yellow area 
denotes the region of inverted charge transport |TMR| > 1. 
The black contours separating two yellow areas, represent the 
divergence of TMR where N p = —TV . Any contour sepa- 
rating blue (red) and yellow represents N° = (N p = 0). 
The black dots indicate the symmetry configurations chosen 
in Fig. H; a), (b) and (c). 



This means that the spin- valve effect in the static linear- 
response regime in lowest order in the tunnel coupling 
cannot be reversed. 

An inverted spin-valve effect was found even with- 
out time-dependent parameters for spin-depedent tun- 
neling^ and for nonlinear bias4^ In contrast, the charge 
transport inversion, which we find for the pumped charge, 
corresponding to |TMR pump | > 1, is not expected in a 
static situation, because it is induced by thermal trans- 
port through the excited state. In the static linear re- 
sponse, thermal transport through an excited state is al- 
ways a minor correction to transport through the ground 
state. In contrast, we showed in this section that due 
to the time dependence of the system parameters one 
can realize certain configurations (tuned by the double 
dot's parameters e and A) in which the contribution from 
transport through the ground state is suppressed in the 
time-averaged pumped charge. 

We finally mention results for the TMR in the weak 
interdot coupling regime, which is not discussed in detail 
here. Also in this case the TMR takes values in the whole 
range [— oo, oo] , which can be shown directly from a study 
of the N-DD-F system in this regime: already if only 
one lead is polarized, we find a reversed spin- valve effect 
as well as reversed charge transport in the presence of 
polarized leads, due to polarization dependent quantum 
charge fluctuations. 

A multiple quantum dot in a carbon nanotube in con- 
tact with ferromagnetic leads was recently realized exper- 
imentally^ - a system which is at the basis of our study in 



Sees. II VI and fVl The application of time-dependent fields 
to such a setup would therefore allow to experimentally 
verify our results regarding spin-dependent transport. 



VI. CONCLUSION 

We have investigated charge and spin pumping 
through a serially coupled double quantum dot in both 
the weak and strong interdot coupling regimes, contacted 
to either normal or ferromagnetic reservoirs. We take 
into account Coulomb interaction with the only restric- 
tion that double occupation of a single dot is excluded. 
When the two dots are weakly coupled, interaction- 
induced quantum charge fluctuations occur, resulting in 
a renormalized level difference. We show that this level 
rcnormalization can be directly measured as a node in the 
pumped charge. In the case of strong interdot coupling, 
we find that hybridized double-dot states contribute to 
the transport with opposite sign. If their level spacing is 
within thermal reach, the time-averaged transport of the 
excited channel can outmatch the ground-state one. This 
results in a sign change in the averaged pumped charge, 
uniquely found in the presence of time-dependent fields. 
Including one ferromagnetic contact, both mechanisms, 
i.e., quantum charge fluctuations in the weak interdot 
coupling and thermal accessibility of the excited state 
in the strong interdot coupling, enable the possibility to 
pump spin in the absence of net charge transport. 
We finally studied the spin-valve effect for non-collinear 
ferromagnetic contacts. Depending on the system's pa- 
rameters, we find a normal as well as an anomalous spin- 
valve effect, and even an actual change in the transport 
direction can be observed, leading to a diverging tunnel- 
ing magnetoresistance. 
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Appendix A: Perturbation expansion for weak 
interdot coupling 

The single-electron eigenstates of the isolated double- 
dot system are given by the bonding and the antibond- 
ing states, see Sec. IIIIBI If the level spacing between 
these two states is large with respect to the coupling to 
the leads T, off-diagonal elements of the reduced density 
matrix become fundamentally important only in second 
order in V (Ref. 1551 ) but are suppressed in lowest order 
in the tunnel coupling, T < k^T . The present work is 
limited to the case of weak coupling, where a description 
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up to first order in T is justified. In Sees. IIII Al and HV A[ 
we are interested in the regime where the splitting of the 
two levels is small and coherent oscillations between the 
related states play an important role already in first order 
in r. Therefore, the interplay of the internal dynamics 
and the dynamics due to the coupling to the leads is most 
interesting when all three parameters are of the same or- 
der, A ~ e ~ r, and thus have to be considered on the 
same level of approximation^ 

In this regime, when coupling between the left and the 
right dot is therefore small, A < T, it is intuitive to de- 
scribe the double-dot system by the states |Ler) and |Rc). 
These are the eigenstates of the fully decoupled single 
dots, and the corresponding basis does not represent an 
exact eigenbasis of the isolated double-dot system. 
The localized states differ from the hybrid bonding and 
antibonding states only when A is taken into account in 
at least first order. On the other hand the Kernel W, 
entering the last term of the generalized master equation, 
Eq. (|5|), describes tunneling to the leads, and therefore 
starts in first order in T. When treating the parameters 
A, e and T on the same footing, for the calculation of 
the Kernel, the eigenstates of the double-dot system are 
therefore taken into account in zcroth order in the small 
parameters and equal |Lct) and |R<r) within the limits 
of the approximation. This additionally justifies the de- 
scription of the dynamics of the system in terms of the 
localized eigenstates, |Lct) and |Rer). 



Appendix B: Static TMR 

We calculate the conductance for the static case [see 
Eq. (|42p ] in the linear-response regime. We focus on 



strong intcrdot coupling A > T and infinite Coulomb 
interactions U and U' in order to compare to the results 
found in Sec. [Vj The linear conductance in the presence 
of polarized leads, G p , is obtained by solving the sta- 
tionary master equation, Eq. (|7a|) for leads with different 
Fermi energies and by calculating the stationary current 
response from Eq. ©. We find for the conductance 



GP =2e 2 /3 J2 f + M 



+ , s r Lt? r Rt? 



?7— b.a 






CT 2 V [/-(£„)]-+ ((A- Pr) B v 

n [f- (tr,)? + B* 



Po , (Bl) 



with 



C = (pi. 



■ Pr) 2 + (plPr) 2 



2 2 

PlPr 



(B2) 



and the static probability of having an empty dot in low- 
est order in tunneling P = (l + 2e~< 3£b +2e _/3ea ) _1 . The 
linear conductance for unpolarized leads G° is found by 
setting the polarizations to zero, making the correction 
term in the parentheses vanish. We insert the linear con- 
ductances found from Eq. (|B1[) into Eq. (|4~2"|) in order 
to find TMR s tat ■ The static tunneling magnetoresistance 
depends on E in contrast to the pumping TMR. Further- 
more, due to the applied bias, the effective magnetic field 
enters in the static case, which it does not in the case of 
pumping for vanishing bias. In contrast to the pumping 
case, TMR s tat takes always values within in the interval 
[0,1]. 
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